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Abstract 

Background: Mpumalanga in South Africa is committed to eliminating malaria by 201 8 and efforts are increasing 
beyond that necessary for malaria control. Differential Equation models may be used to study the incidence and 
spread of disease with an important benefit being the ability to enact exogenous change on the system to predict 
impact without committing any real resources. The model is a deterministic non-linear ordinary differential equation 
representation of the dynamics of the human population. The model is fitted to weekly data of treated cases from 
2002 to 2008, and then validated with data from 2009 to 201 2. Elimination-focused interventions such as the scale-up 
of vector control, mass drug administration, a focused mass screen and treat campaign and foreign source reduction 
are applied to the model to assess their potential impact on transmission. 

Results: Scaling up vector control by 1 0% and 20% resulted in substantial predicted decreases in local infections with 
little impact on imported infections. Mass drug administration is a high impact but short-lived intervention with 
predicted decreases in local infections of less that one infection per year. However, transmission reverted to 
pre-intervention levels within three years. Focused mass screen and treat campaigns at border-entry points are 
predicted to result in a knock-on decrease in local infections through a reduction in the infectious reservoir. This 
knock-on decrease in local infections was also predicted to be achieved through foreign source reduction. Elimination 
was only predicted to be possible under the scenario of zero imported infections in Mpumalanga. 

Conclusions: A constant influx of imported infections show that vector control alone will not be able to eliminate 
local malaria as it is insufficient to interrupt transmission. Both mass interventions have a large and immediate impact. 
Yet in countries with a large migrant population, these interventions may fail due to the reintroduction of parasites 
and their impact may be short-lived. While all strategies (in isolation or combined) contributed to decreasing local 
infections, none was predicted to decrease local infections to zero. The number of imported infections highlights the 
importance of reducing imported infections at source, and a regional approach to malaria elimination. 



Background 

Mpumalanga is one of three malaria-endemic provinces in 
South Africa. As South Africa is committed to eliminating 
malaria by 2018, efforts are increasing in each of these 
three provinces beyond that which was necessary for 
malaria control [1]. In shifting focus from control to 
elimination, additional activities need to be incorporated 
into the operational strategy as elimination may not be 
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achieved through a "more of the same" approach [2]. 
Groups such as the MalERA Consultative Group on 
modelling, have recognized the contribution mathema- 
tical modelling can make to the elimination of malaria 
globally and have highlighted priority areas that modelling 
can inform, such as optimal resource allocation, strategies 
to minimize the evolution of drug and pesticide resistance, 
assessment of new tools to interrupt malaria transmission, 
assessment of combinations of such tools, the coverage 
targets and expected timelines needed to achieve elimi- 
nation goals and the assessment of operational feasibility 
with respect to costs and human resource capacities [3]. 
While such models may be used to understand/explore 
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the underlying system, an important benefit of modelling 
is the ability to enact exogenous change on the system 
to predict impact without committing any real resources. 
This is particularly important in environments with scarce 
competing resources. In this paper, a population-based 
non-linear ordinary differential equation model is used to 
simulate malaria transmission in Mpumalanga province 
to assess the potential impact of various policy interven- 
tions that may be used to achieve malaria elimination. 
Malaria in Mpumalanga has been documented exten- 
sively [4-10]. Sharing borders with both Mozambique and 
Swaziland (Figure 1), Mpumalanga experiences seasonal 
unstable transmission that is prone to sporadic outbreaks 
from the first rains in October to late May. The Ehlanzeni 
District on the eastern border of Mpumalanga is most 
affected by malaria with the number of imported cases 
in Mpumalanga overtaking locally sourced cases in recent 
years (Figure 2). Between 2002 and 2012, 41% of cases 
were sourced in South Africa and 54% sourced from 
Mozambique (the remaining 5% being sourced from other 
African and Asian countries). Source of infection has 
been determined for all cases in the province, whereby 
a case is classified as imported if the patient travelled 
to a malaria-endemic area in the past month or if there 
is no evidence of local transmission (vectors or cases 
within 500 m radius of the place of residence) [11]. The 



proportion of imported cases has increased from 39% 
in 2002 to 87% in 2012 [12]. Extensive vector control 
through indoor residual spraying, the implementation of 
artemisinin-based combination therapy policy of arten- 
sunate plus sulphadoxine-pyremethamine in 2003, fol- 
lowed by artemether lumefantrine (AL) in 2006 and the 
Lubombo Spatial Development Initiative (LSDI) are con- 
sidered responsible for the decline in malaria cases and 
malaria deaths in the province [4]. The LSDI malaria pro- 
gramme was a collaborative project conceived by the 
Medical Research Councils of South Africa, Swaziland 
and Mozambique to decrease malaria in the areas sur- 
rounding the Lubombo Mountains [13]. With well- 
functioning programmes in place in both South Africa 
and Swaziland, intervention took place primarily in 
Maputo province and Gaza Provinces in Mozambique. 
The programme was terminated early in September 2010 
and the resultant reduced IRS in Maputo thereafter coin- 
cides with the increase observed in malaria incidence in 
Maputo from 2011. 

Differential Equation or compartment models, have 
been used in the past to study the incidence and spread 
of disease, and the impact of interventions such as 
drug treatment and parasite control [14]. Compartment 
models and their applications in malaria in particular, 
have a history that spans more than 100 years [15]. 




Figure 1 Mpumalanga province. A map of Mpumalanga Province in relation to Mozambique and Swaziland (Source: Mpumalanga Malaria 
Elimination Programme (unpublished)). 
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Applications of mathematical modelling in Mpumalanga 
include a climate-based fuzzy distribution model of 
malaria transmission in sub-Saharan Africa (including 
a region containing Mpumalanga) [16]. Coleman et al 
used the SaTScan methodology in Mpumalanga to detect 
local malaria clusters to guide the provincial control pro- 
gramme [17], and Montosi et al. considered soil-water 
content as a driver of malaria incidence; applying both an 
ecohydrological model and a transfer function model to 
incidence data in three South African provinces (includ- 
ing Mpumalanga) [18]. The model presented in this paper 
is used to assess the impact of proposed policy interven- 
tions in Mpumalanga. This is the first study designed for 
this purpose in Mpumalanga and the first to do so since 
the call for malaria elimination in South Africa. A deter- 
ministic population-based non-linear ordinary differen- 
tial equation model fitted to the Mpumalanga malaria 
data, is used to predict the impact of the following inter- 
ventions (alone and in combination): scale-up of vector 
control, mass drug administration (MDA), a focused mass 
screen and treat campaign (MSAT) and foreign source 
reduction. 

Methods 

Transmission model 

The model is a deterministic ordinary non-linear dif- 
ferential equation representation of the dynamics of the 
human population. In the model, the population is divided 
into nine compartments: the susceptible population (S), 
the asexual blood stage only (B/ and By ) for locally 
and imported infections respectively and the infectious 
gametocyte stage (1/ and If) for locally and imported 



infections respectively (Figure 3 with parameter descrip- 
tions in Table 1). The blood stage and infectious stage 
compartments are further stratified according to whether 
the infection is treated or not. The liver stage of the 
infection is incorporated as a delay in the flow between 
the susceptible and blood Stage compartments. As this 
is a low transmission environment, immunity and super- 
infection are rare and are excluded from this model. While 
the seasonal nature of transmission is incorporated in the 
model, the mosquito population is not modelled directly 
as it is assumed that the mosquito dynamics operate 
on a faster time-scale than the human dynamics and as 
such the mosquito population can be considered to be at 
equilibrium with respect to changes in the human popu- 
lation [19]. A full description of the model is presented in 
Additional file 1. 

Data fitting 

The model is fitted to weekly incidence data of treated 
cases from 2002 to 2008, and then validated with data 
from 2009 to 2012. Ethical approval for use of the 
data was obtained from the University of Cape Town 
Human Research Ethics Committee and the Mpumalanga 
Department of Health. The seasonal forcing functions 
(functions that determine the seasonal behaviour of trans- 
mission in the area) for local and imported cases are 
derived from the data. Silal et al. describes in detail 
the characteristic triple peaked pattern in the incidence 
data with peaks in the malaria season occurring in 
September/October, December/January and April/May. 
While locally sourced infections exhibit this triple-peaked 
pattern, imported infections occur mainly in the second 
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Figure 3 Model flowchart. Flowchart underlying the population level ordinary non-linear differential equation model of malaria transmission. 
(I - locally sourced infections, f - imported infections, u - untreated infections, tr- treated infections). 



Table 1 Parameter table 



Parameter 


Description 


Value 


Source 


N 


Population size 


4 x 10 6 


[20] 


\x 


Mortality Rate 


105 
10000 


[21] 


8 


Natural recovery period 


26 weeks 


[22-24] 




Period between liver stage and blood stage 


7 days (5-10) 


[25-27] 


G 2 


Period between blood stage and onset of gameto- 
cytemia 


1 week 


[23,28] 


r 


AL elimination half-life 


6 days (3-6) 


[29] 


T 


Time to seek treatment 


1 /2 week 


Expert opinion 


P 


Proportion that receive treatment 


0.95 


[30,31] 


seasi 


Seasonal forcing function for locally sourced cases 


Derived from data 


[12] 


seosf 


Seasonal forcing function for imported cases 


Derived from data 


[12] 


Pi 


Annual number of mosquito bites per person x 
proportion of bites testing positive for sporozoites 


39.170 (38.894,39.448) 


Estimated from model fitting process 




Force of imported infections 


0.002163 (0.002124, 0.002202) 


Estimated from model fitting process 




Force of locally sourced infections 


(1 ^[f])Ax' tu+ ''f^ 




vc[t] 


vecov x vceff 






vecov 


Vector Control Coverage 


0.22-0.90 


Derived from data 


vceff 


Effectiveness of vector control 


0.9060 (0.8884, 0.9212) 


Estimated from model fitting process 



Table providing the values, descriptions and sources of the parameters driving the mathematical model of transmission. 
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two peaks of the malaria season [12]. The two sea- 
sonal forcing functions were derived using the "Seasonal 
decomposition of Time series by LOESS" (STL) methods 
for extracting time series components [32]. In order 
for the data-fitting process to be plausible, interven- 
tions that were implemented between 2002 and 2008 
were included in the model, namely, ACT drug therapy 
and Indoor Residual Spraying (IRS). Ngomane and de 
Jager outline in detail the IRS procedure and physical 
structures sprayed in the province between 2001 and 
2009 [5]. 

The model is run from 1990 to reach a steady state 
before being fitted to data from 2002. The model out- 
put (local and imported treated cases) are fitted to the 
data from 2002 to 2008 using the maximum likelihood 
approach by treating the model output as the rate X of 
the Poisson distribution. The parameters /3/, veff and Ay- 
are estimated through the data fitting process with ini- 
tial values sampled from a Latin square framework. The 
model with the estimated parameter values is then run 
for a further three years to be further validated by com- 
parison to data between 2009 and 2012. The impact of 
routine drug therapy and IRS implemented between 2009 
and 2012 is also included in the model. A full descrip- 
tion of IRS and the data-fitting method are presented in 
Additional file 1. 

Assessing elimination 

Elimination of a disease is a term that has had several 
definitions over time [33]. Currently, the World Health 
Organization(WHO) defines elimination generally to be 
"Interrupting local mosquito-borne malaria transmission 
in a defined geographical area, i.e. zero incidence of 
locally contracted cases, although imported cases will 
continue to occur. Continued intervention measures are 
required" [34]. The framework of deterministic differen- 
tial equation models are such that compartments may 
approach zero but will never actually decrease to zero; 
hence it is technically impossible for the model to pre- 
dict zero incidence of locally contracted cases. It is then 
necessary to set a threshold below which the number 
of locally contracted cases is deemed equivalent to zero. 
This approach has been used in several papers. For exam- 
ple, defined elimination as having been achieved when 
parasite prevalence is reduced to 0.0001% and the rate 
of change in parasite prevalence thereafter is negative 
[22]. Maude et al. defined elimination to be achieved 
when there is fewer than one malaria parasitaemic indi- 
vidual in the population [35]. The threshold for elim- 
ination that is adopted in this paper is less than one 
locally sourced malaria infection per year in the popu- 
lation and the rate of change of locally sourced infec- 
tions is negative thereafter. When the model predicts 
that the number of locally contracted cases is below this 



threshold and the rate of change of locally sourced infec- 
tions is negative thereafter, then elimination is predicted 
to occur. 

Results 

The results of the model fitting and validation are pre- 
sented first before evaluating the predicted impact of the 
elimination-focused interventions. 

Model fitting and validation 

The parameters driving the model and their 95% con- 
fidence intervals estimated through data-fitting pro- 
cedures are presented in Table 1. A parameter that 
is usually unknown and estimated from the data is 
the proportion of cases treated (p). Case data usu- 
ally includes cases that have been treated, compris- 
ing patients presenting themselves at a health facility 
(passive case detection), or those cases that have been 
detected actively. There is often no indication of the 
number of infections that have remained untreated, 
and hence there is no/little data from which to derive 
p. Castillo-Riquelme et al. conducted household sur- 
veys in Mozambique and South Africa between 2001 
and 2002 to evaluate treatment-seeking behaviour for 
malaria-related events [30]. It was found that in the 
Tonga sub-district of Mpumalanga, all of the 457 peo- 
ple with recent cases of malaria (previous month) 
sought treatment, with 98.5% of cases being treated 
at a public health facility. More recently, Hlongwana 
et al. conducted a study on the knowledge and prac- 
tices towards malaria in Bushbuckridge Municipality in 
Mpumalanga in 2008 after South Africa was declared 
ready for malaria elimination [31]. The study revealed 
that 99% of respondents would seek malaria treatment 
(95% Confidence interval: (97.5, 99.5)%) with 82% doing 
so within 24 hours of the onset of symptoms. Based on 
these two studies, a probability of treatment of 95% is 
assumed in the model for the entire modelling period 
as it is informed by studies conducted in 2002 and 
2008. 

Model fitting was performed using different starting 
values of the parameters with the optimization routine 
reaching the global minimum in almost all fits. The nar- 
row confidence intervals of the parameter estimates are 
indicative of the stability of the estimates. Figure 4(a) 
shows the fit of the model to the data for the treated 
cases that were locally and imported while Figure 4(b) 
shows the application of the model to treated case data 
from 2009 to 2012. The model captures the timing of 
the season well. As seen in Figure 4(b), there is a sud- 
den unanticipated rise in the number of imported cases in 
2011 and 2012 coinciding with the end of the Lubombo 
Spatial Development Initiative and the model does not 
capture this. 



Silal etal. Malaria Journal 2014, 13:297 
http://www.malariajournal.eom/content/13/1/297 



Page 6 of 12 



(a) (b) 




2000 2002 2004 2006 2008 2009 2010 2011 2012 2013 



Figure 4 Model fitting and validation to treated case data, (a) Model fitting results to locally sourced (blue) and imported (green) reported 
treated malaria cases and (b) Model validation results to locally sourced (blue) and imported (green) reported treated malaria cases from 2009 to 
2012. 



Interventions 

The model is used to predict the impact of the following 
interventions (alone and in combination): scale-up of vec- 
tor control, mass drug administration (MDA), a focused 
mass screen and treat campaign (MSAT) and foreign 
source reduction. 

Scaling up vector control 

The primary exogenous vector control intervention in 
use currently is IRS with larviciding at specified sites. 
As these interventions act directly on local vectors and 
hence local transmission, their impact is included in the 
model as a percentage decrease in /3/, the number of local 
human contacts with infectious mosquitoes. By 2012, the 
vector control activities decreased $\ to 30% of its initial 
value (Figure 5, black). Scaling up vector control (through 
increased IRS or larviciding for example)so as to decrease 
Pi by a further 10% (Figure 5, red) and 20% (Figure 5, blue) 
results in local infections decreasing substantially with- 
out any sizeable decrease in total infections. This is to be 
expected as the majority of infections are imported, and 
as such vector control may assist with onward local trans- 
mission of these imported infections but not change the 
number of imported infections itself. IRS already takes 
place on a large scale in the Ehlanzeni district [5]. Scal- 
ing up vector control by these amounts will require other 
vector control activities to be implemented such as larvi- 
ciding. To achieve such reductions will require large-scale 



identification and coverage of breeding sites and even still 
the model does not predict the reduction to zero of local 
cases. Local infections will still occur because of onward 
transmission from imported infections. 

Mass interventions 

One intervention that simultaneously impacts both local 
and imported infections is mass drug administration. 
Mass drug administration involves treating all individu- 
als without prior screening to assess disease status. While 
MDA is aimed at the entire population of interest, it is 
rarely the case that every single individual will be treated 
and hence MDA should be modelled with a less-than 
100% coverage. The choice of drug is key to the interven- 
tion as drugs that infer a long period of chemoprophylaxis 
may result in fewer infections after the intervention, but 
also expose parasites to sub-therapeutic levels of drugs 
which may in turn lead to the development of resis- 
tance. The most likely choices for drugs in Mpumalanga 
are the first-line of treatment, artemether-lumefantrine, 
and dihydroartemsinin plus piperaquine with Primaquine, 
inferring a protective period of approximately one month 
[36]. The timing of the MDA is also vital to its effec- 
tiveness. It may be the case that performing MDA at the 
trough of the season will result in fewer malaria infec- 
tions and a decrease in the infectious reservoir, leading 
to fewer infections at the peak of the season. Figure 6 
shows the predicted impact of MDA applied over a two 
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month period annually to the whole population with 
80% coverage. The timing of MDA was investigated at 
both the peak (Figure 6, red) and trough of the sea- 
son (Figure 6, blue). It is predicted that applying MDA 
leads to a substantial decrease in local infections regard- 
less of the timing of application. However, the predic- 
tions show that applying annual rounds of MDA at the 



peak of the season leads to substantial decreases in both 
local and imported infections unlike MDA at the trough 
of the season. A possible reason for this difference in 
impact is the large proportion of imported infections. 
Applying MDA at the seasons trough decreases imported 
infections then but there is no further decrease dur- 
ing the peak of the season as imported infections are 
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Figure 6 Annual rounds of MDA. Predicted impact of annual rounds of MDA performed at the peak (red) and trough (blue) of the season on local 
and imported infections. 
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sourced outside South Africa and so, are not impacted 
by the decrease in onward transmission brought about by 
MDA. What is also evident from Figure 6 is that annual 
rounds of MDA are predicted to be insufficient to bring 
about malaria elimination. A single round of MDA is pre- 
dicted to cause an immediate decrease in infections but 
future transmission recovers quickly. Figure 7(a) shows 
the application of six consecutive two-monthly rounds 
of MDA to the malaria transmission model. In this sce- 
nario, local infections are predicted to decrease substan- 
tially during the MDA, reaching less than 1 local case 
per year (1/52 cases per week) 29 weeks after the start 
of the intervention and take approximately two years to 
recover to pre-MDA levels once the MDA has stopped. 
Imported infections are predicted to decrease substan- 
tially during MDA acting as a form of intermittent preven- 
tative treatment and recover immediately after the MDA 
cycle. 

MDA is a resource-intensive process attempting to 
access an entire population of interest while adminis- 
tering a drug regardless of whether individuals have the 
disease or not. Mass screen and treat campaigns on the 
other hand treat only those that have tested positive 
for the disease. Figure 7(b) shows the predicted impact 
of administering MSAT at the border to residents of 
Mpumalanga who have imported infections i.e. before 
entering Mpumalanga. The rationale behind this inter- 
vention is that it is less resource-intensive than MDA 
and specifically targets imported infections before they 



enter the province and impact local transmission. In 
Figure 7b, MSAT is applied continuously in the model 
to new imported infections before entering Mpumalanga 
for six months from November to April with a coverage 
below 100% as many imported infections may be missed 
for reasons such as illegal immigration and sensitivity of 
the screening tools. Figure 7b also shows the substantial 
predicted decrease in locally sourced infections that can 
be achieved by treating local people who have imported 
infections only through MSAT with 70% coverage. Once 
again however, as soon as the intervention stops, imported 
infections are predicted to revert immediately to previous 
levels while local infections take approximately two years 
to reach previous levels. Lower coverage rates (< 70%) 
have also been explored, with the result of even smaller 
decreases in locally sourced infections. 

Combining interventions 

Figure 8 shows the predicted impact of a number of dif- 
ferent combinations of interventions on locally sourced 
infections. The red line depicts the predicted impact of 
six consecutive two-monthly rounds of MDA followed 
immediately by MSAT (at the border on local people with 
new imported infections) enacted annually for six months 
from November to April at a 70% coverage rate. While 
these two mass interventions are not enough to eliminate 
local infections, a new substantially lower stable cycle is 
reached for local infections. If these interventions are sup- 
plemented with increased vector control (so as to decrease 
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Figure 8 Combination of interventions. Predicted impact of combination of interventions on local infections: Black: No additional interventions, 
Red: 70% coverage of MSAT on local population with new imported infections following 6 months of continuous MDA at 80% coverage, Blue: same 
as red (MDA+MSAT) with increased vector control to decrease transmission by a further 20%, Green: six consecutive two-monthly rounds of MDA 
with increased vector control, and 70% decrease in the foreign force of infection, Purple: six consecutive two-monthly rounds of MDA with zero 
imported infections. 



transmission by a further 20%), local infections reach very 
low levels with a maximum of 16 local infections per week 
at the peak of the season (blue). However, this resource- 
intensive combination of interventions is also predicted 
to be insufficient to achieve elimination of local malaria 
infections. The modelling of MSAT shows that acting 
on the imported infections of local people traveling back 
into Mpumalanga alone decreases the number of locally 
sourced cases. It is expected that the impact would be 
greater if the source of imported infections were targeted. 
Thus the green line is the model prediction based on 
foreign source reduction i.e. if the force of imported infec- 
tion Xj were decreased by 70% and locally six consecutive 
two-monthly rounds of MDA were applied and vector 
control was increased. Additionally, the purple line shows 
the predicted impact on local infections if there were no 
imported infections at all. In this scenario, the only other 
intervention that was applied was six consecutive two- 
monthly rounds of MDA. This combination is sufficient to 
eliminate local malaria in Mpumalanga. It is predicted to 
take only 28 weeks on average from the end of the MDA 
cycle to reach below one local infection per year. 

Discussion 

South Africa has been employing vector control to con- 
trol malaria since 1931 [37]. Reliance on vector control 
has been such that insecticide resistance to pyrethroids in 
2000 resulted in a surge in malaria cases that could not 
be controlled through drug therapy alone. Consistent and 
large-scale IRS is considered one of the key reasons why 



that malaria has been so well controlled in the country 
and in Mpumalanga. Scaling up IRS even further through 
targeted larviciding of vector breeding sites may be con- 
templated as a strategy to achieve malaria elimination. 
Scaling up vector control in the mathematical model for 
Mpumalanga, allowed local infections to decrease to a 
new equilibrium but did not eliminate local infections 
or decrease total infections substantially. It is expected 
that total infections will not decrease substantially as 
the majority of infections in Mpumalanga are imported 
and local vector control will not impact the number of 
imported infections in the province but will work towards 
decreasing onward transmission of these infections. It is 
because of this constant influx of imported infections that 
vector control alone will not be able to eliminate local 
malaria. A "more of the same" approach does not appear 
to work as it is insufficient to interrupt the transmission 
stable cycle. 

Mass drug administration is a resource intensive inter- 
vention that needs to be acted out quickly, systemati- 
cally and efficiently. It is a strategy that has not found 
favour in recent years; one reason being that drug resis- 
tance is a feared consequence of MDA [38]. In a global 
environment seeking to protect the artemisinins (and 
their partner drugs) from drug resistance, the choice of 
drug for mass interventions has both global and eco- 
nomic importance. Economically, the drug needs to be 
affordable to be deployed on a large scale and glob- 
ally, if resistance spreads to a point where the drug is 
no longer useful, the drug should be one that can be 
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sacrificed. While the emergence of drug resistance has 
not been directly linked with MDA, it is possible that 
MDA will increase the selection pressure on the para- 
site population, with the possibility of losing the drug 
eventually [39]. In contrast, in fitting a mathematical 
model to trial results from Western Cambodia to assess 
the effects of elimination strategies and their interactions 
with artemisinin resistance, Maude et al. found that the 
proportion of artemisinin resistant infections increased 
quickly when ACT was introduced for treatment (in an 
area where treatment comprised non-artemisinin anti- 
malarial drugs and a low level of artemisinin monother- 
apy) and only slightly more when MDA with an ACT was 
implemented [35]. 

Both MDA and MSAT are interventions that have 
an immediate and large impact. Yet in countries with 
a large migrant population or populations with sub- 
optimal coverage, these interventions can fail because 
of the reintroduction of parasites and their impact 
may be short-lived [39]. Given the large proportion 
of imported cases in Mpumalanga, the model predicts 
that when a mass intervention targeting imported infec- 
tions is stopped, both local and imported infection lev- 
els revert back to pre-intervention levels with a few 
years. This is because forces within Mpumalanga do 
not determine the level of imported infections; it is 
determined by the prevalence of malaria in the source 
country itself. Hence as the results of applications of 
MDA and MSAT in the model have shown, the large 
predicted impact can only be sustained if the mass 
interventions are applied often, and this is a resource- 
intensive strategy. Applying MSAT to local residents with 
imported infections in the model has predicted a sub- 
stantial "knock on" decrease in local infections. This 
is expected as decreasing imported infections decreases 
the infectious reservoir in the province, which in turn 
decreases local malaria transmission. The model also 
predicted the substantial impact decreases in the for- 
eign force of infection (foreign source reduction) has 
on local transmission. The scenario of zero imported 
infections was the only situation in which the model 
predicted that elimination of malaria could be achieved. 
This highlights the importance of source-reduction, mon- 
itoring imported infections and the receptivity of key 
areas within Mpumalanga. If the malaria vectors are 
present and ecological and climatic factors favour trans- 
mission in these areas, then onward transmission is 
probable even in the presence of good malaria con- 
trol [40]. One strategy to decrease malaria transmis- 
sion to zero would be to eliminate the mosquito vec- 
tor population. This is highly unlikely to be feasible. 
As a result, even if low levels of malaria prevalence 
have been achieved, imported infections will augment 
the infectious reservoir, and since the vector remains, 



imported infections may lead to onward transmission 
to the local population and a resurgence of malaria 
generally [2]. 

In this deterministic model, a threshold was defined to 
determine elimination as it is mathematically impossible 
to achieve zero cases in a differential equation frame- 
work. As is the case with selecting any threshold, more 
severe thresholds could always be selected. It is of inter- 
est in this manuscript that malaria elimination could 
not even be predicted with the current threshold, let 
alone more severe thresholds. In a review on the his- 
torical and current definitions of malaria, Cohen et al 
defines three states of malaria transmission to be con- 
trolled low-endemic malaria, elimination and controlled 
non-endemic malaria. This third state describes the situ- 
ation where the interruption of endemic transmission has 
occurred but there is still malaria resulting from onward 
transmission from imported infections and this onward 
transmission is sufficiently high that elimination has not 
yet been achieved [33]. This implies that if all onward 
infection from imported infections could be prevented, 
elimination of malaria would follow naturally. The results 
of this modelling exercise suggest that Mpumalanga is in 
this third state of transmission as elimination of malaria 
is only predicted to be possible with unrealistic, resource- 
intensive interventions that result in a drastic reduction in 
imported infections. In realizing that the key to decreasing 
local infections further is to prevent imported infections, 
new approaches must be explored both nationally and 
regionally. 

The transmission model presented in this paper is a 
deterministic population-level one considering only the 
population of Mpumalanga. While population-level mod- 
els are useful to assess aggregate effects, the entire pop- 
ulation is treated as homogenously affected by malaria 
transmission. Thus there is no scope to include spatial 
variation or heterogeneous behaviour in the transmission 
model or in the interventions themselves. Imported infec- 
tions play an important role in transmission dynamics in 
Mpumalanga and human migration is modelled indirectly. 
Current work includes disaggregating the population 
into smaller groups such as administrative districts, and 
using a stochastic, meta-population model of transmis- 
sion to model these districts while explicitly incorporat- 
ing human migration between Mozambique (the main 
source of imported infections) and Mpumalanga. This 
will allow one to assess the impact of deploying interven- 
tions sequentially in the administrative districts as well 
as in Mozambique itself. For example, Silal et al found 
that infections in Bushbuckridge municipality occur at 
higher levels at the start of the season compared to other 
municipalities in Mpumalanga [12]. It may be of inter- 
est then to assess if intervening in this municipality at 
the start of the season has a knock-on effect on malaria 
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in the rest of the province. Spatial heterogeneity will also 
allow for improved modelling of interventions that tar- 
get imported infections such as border control through 
MSAT, and source reduction in Mozambique. Future work 
also includes incorporating vector population dynamics 
in the model so that vector control activities such as 
indoor residual spraying and larviciding may be modelled 
explicitly. 

Conclusions 

South Africa aims to achieve malaria elimination by 2018. 
This requires synchonized action in the three provinces 
in which malaria occurs in order to decrease local infec- 
tions to zero. In the case of Mpumalanga, given the large 
proportion of imported cases, interventions also need 
to target imported infections to decrease the infectious 
reservoir impacting local transmission. In isolated coun- 
tries, a nationally focused elimination programme may 
stand a better chance of success than countries with high 
levels of visitation from higher transmission regions [41]. 
This paper has used population level mathematical mod- 
elling to model transmission in Mpumalanga and test out 
strategies (MSAT, MDA, increased vector control and for- 
eign source reduction) that may be used to achieve malaria 
elimination. While all strategies (in isolation or com- 
bined) contributed to decreasing local infections, none 
was able to decrease local infections to zero due mainly to 
the continuous stream of imported infections highlight- 
ing the importance of source reduction and a regional 
approach to malaria control and elimination. Disaggrega- 
ting the model into smaller groups will allow for the spatial 
heterogeneity required to optimize elimination strategy 
that may lead to different results. Mathematical modelling 
has the potential to inform government policy to achieve 
malaria elimination and with effective and efficient inter- 
ventions, adequate sustainable finance, local and inter- 
national political commitment and an epidemiological 
understanding of malaria elimination, malaria elimina- 
tion in Mpumalanga may be possible in the foreseeable 
future. 
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